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ON  THE  UTILITY  OF  THE  m-6-8  POTENTIAL  FUNCTION 


t  tt 

H.    J.    M.    Hanley     and  Max  Klein 


The  calculation  of  transport  and  equilibrium  properties  of  simple  gases  with  the 
m-6-8  potential  function  is  discussed.     Properties  referred  to  specifically  are  the 
viscosity,    diffusion,    thermal  conductivity  and  second  virial  coefficients,   and  the  thermal 
diffusion  factor ,     Gases  referred  to  are  argon,    krypton,   xenon,    nitrogen,    oxygen,    carbon 
dioxide  and  methane.     We  also  discuss   in  detail  the  genera]  behavior  of  a  model  potential 
with  respect  to  the  fitting  of  data  and  the  selection  of  potential  parameters.      The  m-6-8 
appears  to  be  the  simplest  model  potential  which  can  satisfactorily  be  used  to  fit  data  over 
a  wide  temperature  range.     The  relationship  of  the  potential  to  independent  theory  is 
briefly  examined. 

Key  words:     Dispersion  coefficients;   m-6-8  potential;  parameter   selection;   second 
virial  coefficients;  transport  properties. 
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1.     Model  Potential  Functions 

In  this  paper  we  discuss  a  new  model  potential  function,    the  m-6-8,   and  demon- 
strate how  it  can  be  used  to  correlate  and  predict  experimental  data.     Before  going  into 
any  details,   however,   we  would  like  to  review  very  briefly  the  general  development  of 
model  functions  and  their  relationship  to  experiment  and  statistical  mechanics. 

Roughly  speaking,    the  evolution  of  model  potential  functions  parallels  theoretical 
developments   in  molecular   structure  and  inter  molecular  forces.      For  example,    the 
simplest  concept  of  molecular   structure  is  that  of  the  hard  spherical  molecule  and  the 
corresponding  potential  accounts  only  for  repulsion  between  molecules  on  contact.     The 
potential  has  one  parameter,    the  effective  diameter  of  the  spheres.     While  this  concept  is 
obviously  too  simple  --it  cannot,   for  example,    explain  the  stability  of  solids  and  liquids, 
or  give  a  temperature  dependence  to  the  second  virial  coefficient  --    it  was  nevertheless 
very  useful,   and  remains  so  today.     The  hard  sphere  model  can  be  modified  to  incor- 
porate a  simple  inverse  power  repulsive  term  consistent  with  the  picture  of  a  "soft" 
repulsive  spherical  molecule  whose  diameter  decreases  as  the  energy  of  the  molecular 
system  increases. 

Realistic  potentials  have  to  consider  both  the  repulsion  and  attraction  between 
molecules.     This  could  be  taken  into  account  by  using  inverse  power  potentials  of  the  type 

$(r)     =    A/r       -   C/r  (1) 

generally  associated  with  the  name  of  Lennard- Jones.     Here     $  (r)     is  the  potential,   a 
function  of  inter  molecular   separation,     r,     while     m    and    n,      (m.>n)  are  fixed  at  arbitrary 
values,   and    A     and     C    are  parameters  to  be  found  from  experiment.     One  method  for 
finding  such  parameters  is  described  below. 

Several  authors  have  attempted  to  be  more  specific  by  predicting  values  for     m 
and    n.     For  example,   details  of  the  attraction  between  m.olecules  were  developed  fron:i 
quantum  mechanical  perturbation  theory,    i.e.  ,   attraction  was  seen  to  be  equivalent  to  the 
concept  of  induced  molecular  polarizability.     Specifically,   quantum  mechanics   shows  that 
the  leading  term  of  the  induction  effect  varies  as  the  inverse  sixth  power  of    r     thus 
suggesting  that  the  attractive  index,     n    of  eq  (1),    should  be  taken  as  6.     It  is  not  possible, 
however,    to  come  up  with  a  satisfactory  a  priori  value  for  the  repulsive  index    m    from 
quantum  mechanics,    so     m     is  often  taken  to  be  12  for  mathematical  convenience.     With 
this  choice  of    m,     we  have  the  12-6  or,   as  it  is  usually  known,    the  Lennard- Jones 
potential: 

I  (r)     =    A/r        -   C, /r  (la) 


where    A    and     C,    --   the  coefficients  of  repulsion  and  inverse  sixth  attraction, 
respectively  --  are  two  parameters  to  be  found  fronn  experiment  as  before.     This  potential 
has  been  the  standard  for  many  years  and  still  serves  as  the  model  potential  for   many 
molecular  theories. 

In  the  last  fifteen  years  or   so,    several  three-parameter  potential  families  have 

been  proposed  which  try  to  add  flexibility  to  the   12-6.     Most  often,    the  inverse  sixth  term 

is  retained  for  the  attraction  but  an  extra  parameter  is  introduced  for  the  repulsion.      The 

simplest  example  is  the  m-6  family  of  potentials  with  the  form  of  eq  (la)  but    m    replacing 

the  index  12  as  a  parameter;     m    has  then  to  be  found  from  experinnent  along  with    A     and 

C,  •     Other  examples  of  such  three  parameter  potentials  are  the  Kihara,    the  exp:6,   and 
6 

the  Morse  potentials  [l]  .  Unfortunately,  the  application  of  the  three- parameter  fannilies 
to  experiment  often  gave  misleading  or  even  erroneous  results:  thus,  there  seemed  to  be 
no  definite  criteria  to  judge  how  effectively  a  potential  represented  a  given  property  --in 
naany  cases,  even,  it  was  not  clear  if  the  introduction  of  the  extra  parameter  led  to  a  fit 
between  theory  and  experiment  which  was  an  improvemient  over  the  corresponding  fit  with 
the  12-6.  Consequently,  it  also  was  difficult  to  say  which  of  the  new  three-paranneter 
potential  families  was  the  most  generally  useful. 

To  try  to  clear  up  some  of  the  confusion,   we  undertook  a   study  [1-3]  to  investigate 
the  relationship  between  certain  experimental  properties  and  the  potential  function  as 
required  by  statistical  mechanics.     The  nnain  results  of  the   study  are  of  interest  to  the 
present  work  so  they  are  discussed  briefly  here. 

First  of  all,   we  found  that  temperature  ranges  exist  for  different  properties  in 
which  these  properties  are  only  sufficiently  sensitive  to  the  details  of  the  function  to 
produce  information  about  two  parameters:     in  other  words,    in  these  temperature  ranges, 
the  properties  are  insensitive  to  any  extra  parameters  which  might  be  introduced.     For 
equilibrium  properties  this  range  is,    in  reduced  form,   approxinnately 

2.t;^_^.10  (2a) 

and  for  transport  properties  the  range  is  approximately 

2  .  T*^_^    .  5  (2b) 


Numbers  in  brackets  refer  to  references  given  at  the  end  of  the  paper, 
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where    T         ,   -  T / {e         ,/k)     with    T     the  temperature  in  kelvin,     k    Boltzmann's  constant 
and     e^^    /     the  maximum  potential  energy  of  attraction  for  the  Lennard- Jones,    12-6 
potential  [4] . 

The  definition  of  the  insensitive  temperature  range  enabled  us  to  give  a  quantitative 
significance  to  the  terms  "high"  and  "low"  temperatures  --  terms  which  are  important  in 
discussions  of  the  limiting  behavior  of  potential  functions.     Thus,   a  high  temperature  is  a 
temperature  above     T  ?-  10    for  equilibrium  properties  (or  5  for  transport  properties) 

while  a  low  temperature  is  a  temperature  below     T         ,  r^  2. 

Note  that  the  insensitive  temperature  range  covers  the  normal  experimental  range 
for  most  simple  substances   (e.g.  ,    the  second  virial  coefficients  of  argon,    oxygen,    or 
nitrogen  are  insensitive  to  more  than  two  parameters   in  the  potential  function  in  the  tem- 
perature range  200  <  T(K)  <  1000).     In  our  opinion,    many  authors  were  not  familiar  with 
the  insensitive  temperature  ranges  and  this  lack  of  awareness  led  to  most  of  the  confusion 
in  dealing  with  three-parameter  functions.     For  example,   use  of  data  in  the  insensitive 
range  alone  often  forces  a  decision  in  favor  of  one  potential  function  based,    in  the  final 
analysis,    on  the  particular  methods  used  for  treating  the  experimental  errors.     The 
selection  of  the  potential  function  can  therefore  be  made  on  a  totally  irrelevant  basis. 

As  a  second  result  of  our  study,   we  found  that  all  three-parameter  function 
families  are  essentially  equivalent  as  far  as  fitting  low  density  data  is  concerned.     This 
particular  result  can  be  used  to  advantage  since  one  can  choose  any  convenient  three- 
parameter  family,    say  the  m- 6 ,    to  be  typical  of  all  three-parameter  functions.     We  have 
often  done  this  in  our  previous  work  [3]. 

A  third,   and  more  negative,    conclusion  drawn  was  the  clear  recognition  that  the 
three-parameter  functions,   while  an  improvement  over  the  two- parameter   12-6,   were 
still  not  flexible  enough  to  fit  data  as  one  would  like.     In  particular,    they  failed  to  satisfy 
two  important  requirements: 

(1)  The  ability  to  correlate  data  for  a  given  substance  and  property 
at  high  and  at  low  temperatures  with  the  same  set  of  parameters. 

and 

(2)  The  ability  to  correlate,    for  a  given  substance,    two  different 
properties  using  the   same  set  of  parameters  in  a  sensitive 
temperature  range. 


Other  multiparameter  functions,    where  the  parameters  are  related  to  theory  as  far 
as  possible,    have  been  proposed  by  several  authors.      One  recent  attennpt  is  that  of 
Barker  and  Pompe  [5].     Essentially,    their  potential  is  based  on  a   series  of  exponential 
ternns  for  the  repulsion,    which  are  determined  empirically,   along  with  a  series  of  inverse 
powers  for  the  attraction,    which  are  estimated  from  quantum  mechanical  dispersion  theory. 

An  alternative  technique  for  developing  a   multiparameter  function  is  to  start  with  a 
large  drawing  of  a  reasonable  potential  and  then,   using  tabular  values  read  off  the  drawing, 
fit  different  properties  simultaneously.      By  modifying  the  potential  graphically  in  some 
consistent  manner  until  the  best  overall  agreement  between  theory  and  experiment  is 
obtained  [6,7],   a  function  might  possibly  be   selected  out  whose  shape  closely  approaches 
that  of  the  "actual"  potential  function.      See  in  particular,    the  work  of  Dymond  and 
Alder  [7]. 

This  covers  a  brief  outline  of  the  potential  models  used  up  to  the  time  this  work 
was  undertaken.     Very  recently,   potential  functions  have  been  determined  from  low  energy 
scattering  experiments  and  from  spectroscopy.     See,   for  example,    reference  [99]. 

2.     The  Four- Parameter  ,    m-6-8.    Potential 

Although  generally  successful,    the  complex  potentials,    such  as  the  multiparameter 
semi- theoretical  potential  of  Barker  and  Pompe,    or  the  numerical  potential  of  Dymond 
and  Alder,   are  very  specialized.      These  potentials  are  derived  for  a  particular   substance 
and  are  not  always  easy  to  generalize  to  other   substances.      It  is  also  not  easy  to  take  into 
account  possible  changes  either  in  the  assunnptions  or  in  the  data  on  which  these  functions 
are  based.     We  felt  that  one  had  to  connpromise  between  these  sonnewhat  inflexible  complex 
potentials  and  procedures  on  the  one  hand,   and  the  more  flexible  but  over-simplified 
three-parameter  functions  on  the  other  hand.     Our  approach  is  straightforward:     we  con- 
sider the  addition  of  only  one  parameter  to  the  three-parameter   models  used  previously 
and,   furthermore,    select  a  paranneter  with  a  reasonable  physical  basis. 

Extension  of  quantum  mechanical  perturbation  theory  indicates  that,   for   separated 
molecules,    the  attractive  part  of  the  inter  molecular  potential  function  should  be  described 
by  a  series  of  inverse  powers  of  the  intermolecular  distance,    the  first  term  of  which 
depends  on  the  inverse  sixth  power  of  the  intermolecular   separation,   as  already  men- 
tioned.    The  next  term  of  the  series  depends  on  the  inverse  eighth  power  and  is,    therefore, 
a  logical  choice  for  an  additional  term  for  a  model  potential.      Theoretical  estimates  have 
been  made  of  the  relative  magnitudes  of  the  coefficients  of  the  inverse  eighth  and  inverse 
sixth  power  terms  for  a  nunnber  of  substances  [107],    which  enables  one  to  introduce  an 
inverse  eighth  power  dependence  in  terms  of  the  coefficient  of  the  inverse  sixth  power 
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term  without  additional  parameterization.      However,    since  the  validity  of  these  theoretical 
estinnates  has  not  been  clearly  established,    this  cannot  be  done  with  certainty.     It  seems 
reasonable,    therefore,    to  add  an  inverse  eighth  power  term  but  to  take  the  corresponding 
coefficient  as  an  additional  potential  parameter,    thereby  defining  a  four-paranaeter 
function  fanaily. 

The  particular  choice  of  a  three-parameter  potential  function  on  which  to  add  an 
inverse  eighth  power  attraction  is  arbitrary  because,   as  we  have  pointed  out,   all  of  the 
three-paranrieter  functions  are  equivalent  from  the  point  of  view  of  correlating  data.     The 
equivalence  shows  up  especially  in  the  arbitrary  choice  of  the  repulsive  terms.     One  could 
argue,   however,   based  on  our  knowledge  of  high  energy  scattering  data,    that  a  shielded 
coulonnb  form  for  this  repulsion  would  be  most  suitable.     Nevertheless,    we  decided  to 
continue  to  use  the  single  term  inverse  power  repulsion  contained  in  the  m-6  function. 
This  choice  was  made  because   (a)   such  a  form  is   sinaple,   and   (b)  adding  an  inverse  eighth 
power  to  an  m-6  potential  is  a  natural  extension  of  the  work  we  already  have  done  with  the 
nn.-6 . 

2.  1     Definition  of  the  Reduced  m-6-8  Potential 

The  definition  of  the  m-6-8  potential  follows  fronn  eq  (1): 

l(r)     =    A/r""-   C^/r^  -   Cg/r^  (3) 

where     m    is  the  repulsive  index,   and    A,     C ,  ,     and     C       are  the  coefficients  for  repulsion 

D  o 

and  for  the  sixth  and  eighth  power  attraction  terms  respectively.     We  introduce  the  more 
convenient  reduced  form  of  potential  (3)  as  follows. 

If    r         denotes  the  intermolecular   separation  at  the  maximum  attractive  energy, 
m 

then  we  define     $  (r      )     =     -e.     Using  that  at    r      ,     d$/dr     =     0,     one  can  show  that 
m  m 
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with    r'  -  r/r      .     We  now  define  a  parameter,     y,     characteristic  of  the  inverse  eighth 
term, 

Y    =     -Cg/er^  (7) 

thus 

The  potential  is  given  in  eq  (8)  in  terms  of    r',     i.e.  ,    scaled  by    r      ,     but  it  is  frequently- 
more  convenient  to  work  in  terms  of  the  effective  hard  sphere  diameter,     a,     defined  by 
$  (r)     -     0.     Introduction  of  the  definitions 


d    =     r      /a  and  r*    =     r/o  (9) 


m 


then  converts  eq  (8)   into 

,,  (,.,)   __  Ii_L|ll  (A.y  _  l^-yi^-^)]  (^f  _    (Af  (10) 

m-6         \r*/  m-6  Vrv/  \r*/ 

where  we  have  defined     f^l^  (r=;=)     as     l(r*)/e.     Potentials   (8)  and   (10)  are,    of  course, 
equivalent.      Care  must  be  taken,    however,   when  converting  from  one  form  to  the  other  to 
ensure  that  the  factor     d     is  used  properly.     We  will  work  with  potential  (10). 

The  m-6-8  potential  is  sketched  in  figure  1.  Figure  la  shows  how  the  attractive 
part  is  affected  by  an  increase  in  m  with  y  =  O-O-  Figure  lb  shows  how  the  attractive 
part  is  affected  by  increasing  y  with  m  constant.  Note  that  increasing  y  at  constant 
m  has  roughly  the  same  effect  as  increasing  m  at  constant  y.  [We  have  not  shown  the 
associated  small  changes  in  the  repulsive  part  of  the  potential,  such  changes  being  very 
small  when  compared  to  the  changes  in  the  attractive  part.] 

3.      Testing  The  m-6-8  Potential 

We  next  outline  our  approach  for  testing  the  potential.     We  clarify  what  we  propose 
to  do:    Any  potential  function  obtained  fromi  fitting  data  alone  is  only  the  "best"  function 
associated  with  a  particular  functional  form,   a  particular   set  of  experimental  data,   and  a 
particular  criterion  of  best  fit.     A   successful  potential  function,    from  a  correlation  point 
of  view,    is  therefore  not  necessarily  a  valid  one  from  a  theoretical  viewpoint,   although 
the  reverse  must  always  be  true,    i.  e.  ,    a  valid  theoretical  potential  must,    of  necessity, 
correlate  good  data.     When  testing  with  experimental  data,   we  therefore  refer  to  testing 


the  utility  of  a  function  rather  than  testing   its  validity.     Of  course,    if  the  potential 
parameters  are  then  found  to  be  compatible  with  independent  theoretical  estimates,    the 
potential  takes  on  a   significance  beyond  its  ability  to  fit  the  data. 

Substances:     Since  the  potential  function  assumied  is   spherical,    the  necessary  first 
test  has  to  be  tnade  with  nnonatomic  molecules  because  the  potential  function  can,    in 
principle,   describe  the  interactions  between  such  molecules  without  further  approximation. 
[Non- sphericity  not  only  involves  the  potential  function,    but  also  a  description  of  the 
interactions  between  the  internal  degrees  of  freedom  of  the  molecules.]     We  exclude 
quantuni  fluids  for  now.     Hence,    initial  tests   should  be  made  with  argon,    krypton  and 
xenon.     Of  these,    argon  is  the  most  convenient  because  argon  data  for   several  properties 
are  available  outside  the  reduced  insensitive  temperature  ranges  defined  by  eq  (2).     Since 
this  is  not  always  the  case  for  the  heavier  rare  gases,    much  of  this  paper  will  refer  to 
comparisons  between  predicted  and  experimental  values  for  argon. 

Properties:     There  are  a  number  of  fundamental  theoretical  expressions  which 
give  macroscopically   measurable  properties   in  terms   of  potential  functions  and  which  can, 
in  principle,    be  used  to  test  the  utility  of  such  potentials.      Only  a  few  of  these  are 
practical,    however.      The  reason  is   simply  that  the  corresponding   experimental  data   must 
be  available  and,   furthermore,    must  be  precise  enough  to  allow  for  a  choice  to  be  made 
annong  potential  functions  [1-3,9].     In  practice,    the  usable  expressions  are  those  for  the 
viscosity  coefficient,     Tj,     and  the  second  virial  coefficient,     B.     Of  the  two,    the  viscosity 
is  usually  to  be  preferred  since  the  insensitive  range,    eq  (2),    covers  a  smaller  range  of 
temperature  than  does  the  corresponding  range  for  the  second  virial.     Thus,   data  for  a 
given  substance  covering  a  given  experimental  temperature  range  are  nnore  likely  to 
overlap  both  high  and  low  temperature  regions  for  the  viscosity  coefficient  than  for  the 
second  virial  coefficient.     This  must,    however,   be  balanced  against  the  fact  that,    of  the 
two  properties,    the  second  virial  coefficient  is  the  more  sensitive  to  the  potential  at  low 
temiperatures. 

The  expressions  connecting  these  quantities  with  the  potential  function  are  [lO]: 
Viscosity  (to  the  first  Chapman- Enskog  approximation  [ll]) 


'I  16  2      (2,2)* 

rra     Q  (T-'-) 


Second  Virial  Coefficient 

B    =     b     B*  (T*) 
o 

where  (IZ) 

2  3 

b     =    —    TT  Na 
o         3 

(2    2)* 
In  eqs   (11)  and  (12),    the  intermolecular  potential  function  enters  via     Q     '  (T*)     --     the 

collision  integral  accounting  for  the  dynamics  of  a  binary  collision  in  the  viscosity 

process,   and     B*  (T*)     --     the  reduced  second  virial  coefficient.     For  a  given  potential 

(2,2)- 
function,      §(r),    both     Q  (T*)    and     B*  (T*)     can  be  calculated  as  functions  of    T*, 

where     T*  =  T/(e/k)     with     e     the  maximum  well-depth  for  the  given  potential.     Also,    in 

these  equations,     m     is  the  mass  of  a  naolecule,     N     is  Avagadros'   number  and     a    has 

been  defined  previously. 

Data:     The  argon  data  selected  for  test  purposes  [12]  were  extracted  from  the 
following   sources:    viscosity  [13-16];   second  virial  [17].     We  have  carefully  evaluated 
these  data  and  they  appear  to  be  the  best  available  at  this  time.      There  are  problems, 
nevertheless:    values  of  both  the  viscosity  and  second  virial  coefficients  at  the  lower  tenn- 
peratures     (T  <  100  K)  remain  uncertain.     Also,    the  reliability  of  the  published  viscosity 
coefficients  at  high  temperatures  has  been  questioned  recently.     [This  latter  problem, 
however,    now  seems  resolved  [12].] 

Collision  Integrals  and  Reduced  Second  Virial  Coefficients:     Collision  integrals 
and  reduced  second  virial  coefficients  have  been  computed  for  the  m-6-8  potential  with 
9  s  m  <  18  for  values  of    y     such  that    0  <  [m  -  y('^  -   8)].     The  calculation  procedure  and 
tables  are  available  in  another  publication  [8]. 

3.  1     Fit  to  a  Single  Property:     Satisfaction  of  Requirement  (1),    Section  1 

We  began  our  investigation  by  studying  how  the  m-6-8  potential  can  be  used  to  fit 
a  single  property,    namely  the  viscosity.      The  dilute  gas  viscosity  coefficients  were 
examined  separately  at  low  temperatures,   at  high  temperatures,   and  at  all  temperatures. 

Several  techniques  are  available  for  optimizing  potential  function  parameters. 
For  example,    one  can  select  the  parameters  which  give  the  best  fit  of  theory  to  experime 
according  to  least  squares.     We  use  a  graphical  technique,    however,    which  we  have  four 
to  be  both  successful  and  instructive.     To  illustrate  this  approach,   let  us  first  consider 
argon  viscosity  data  in  a  given  temperature  range. 
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In  our  fitting  procedure,   a  potential  with  particular  values  of    m    and    y     is  first 
selected,    i.  e.  ,    only     e     and     a    are  to  be  varied   initially,      a     is  then  also  fixed   initially  at 
some   sensible  preliminary  value.      This  leaves  only     e     to  be  determined.      Let  us  consider 
a  particular  experimental  temperature     T.      Given  exper innental  data  for     T\    and     T     and  the 

initial  choice  of     a,      eq   (11)   can  be  used  to  obtain  an  experimental  reduced  collision  inte- 

(2    Z)* 
gral     Q      '  at  that  value  of     T.      Our  purpose   i.s   to   make  a  correspondence  between  such 

"experimental  collision  integrals"   obtained  at  the  exper  innental  tennperature  s  and  calcu- 
lated from  (11),   and  those  calculated  theoretically  (as  in  [8])  for  the  particular  choice  of 
m    and    y    at  theoretical  temperatures     T*.     The  correspondence  results   in  a  value  of     e /k 
for  the  given  experimental  temperature     T    from  the  relation  between    T    and     T*.     This  is 

accomplished  by  interpolating   in  the  calculated  set  listed  as  a  function  of    T*    for  the     mi 
and    Y    values  chosen  for  the  value  of  the  "experimental  collision  integral"  found  using   (11). 
Using  the  value  of    T*     corresponding  to  this  value  of  the  theoretical  collision  integrals, 
the  experimental  temperature     T    and  the  relation     e/k  =  T/T*,   a  value  of     s/k    is  obtained 
for     T.     Proceeding  in  this  way  for  all  of  the  experimiental  temperatures,   we  then  obtain 
e/k    as  a  function  of    T.     If,    now,    the  value  taken  for     a    is  changed,   a   second  set  of  values 
of     e/k    is  obtained  as  a  function  of    T.     Proceeding  in  this  way  a  set  of  curves  is  obtained, 
one  for  each  value  of    a.     These  curves  represent  the  variation  of     e/k    with    T    and    a 
for  the  given  choices  of     m    and     y.     Now,    if  the  potential  function  were  valid  and  the 
"correct"  value  of    a    used,      e/k    would  not  vary  with  temperature.     Hence,   a  reasonable 
criterion  for   selecting  the  best    a    value  is  to  choose  that  one  which  is  associated  with  the 
least  variation  of     e/k    with    T.     The  appropriate  value  of     e/k    to  take  is  one  within  its 
range  of  variation  with    T    for  that    a.     Thus,    this  procedure  leads  to  a  best  pair,      e/k 
and     a,     for  the  particular     m    and    y.      By  repeating   the  process  for  different     vn,     y 
pairs,    one  can  select  the     a,    e/k,     m,     y     combination  which  gives  the  smallest  overall 
variation  of     e/k    with  temperature. 

Results 

Both  at  high  and  at  low  temperatures,   fitted  separately,    we  were  able  to  find 
several  sets  of     m,     y    and     a    which  gave  a  constant     e/k     (constant  within  experimental 
error).      This   is  just  the  kind  of  result  found  with  the  three- parameter  functions.      We 
were,   however,   able  to  reduce  the  possibilities  to  a  single  set  of     m,     y>     ^    ^Y  requiring 
that  the  same  values  of    a    and     e/k    be  applicable  to  both  high  and  low  temperatures. 

Values  thus  obtained  for  argon  were     m=ll,     y  =  3'0'     a=3.29A,      e/k=153K 
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(1   A  =  10  m). 

This  result  was  encouraging   since  it  meant  that  the  m-6-8  potential  can  satisfy 
requirement  (1)  of  Section  1   --  the   simplest  function  found  thus  far  to  do  so.      It  is 
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important  to  remark,   however,    that,    in  order  to  do  this,    data  are  required  for  both  the 
high  and  low  temperature  ranges.      In  other  words,   a  fit  of  data  at  low  temperatures  only 
is  not  necessarily  applicable  to  data  at  high  temperatures,   and  vice  versa. 

3.2     Simultaneous  Fit  of  the  Viscosity  and  Second  Virial  Coefficients: 
Satisfaction  of  Requirement  (2),    Section  1 

In  a  previous  paper  [2]  we  mentioned  another  graphical  and  interpolation  scheme 
whereby  a  simultaneous  fit  of  the  viscosity  and  second  virial  coefficient  for  a  given  gas 
could  be  used  to  give  potential  parameters  and  which  would,    at  the  same  time,   demon- 
strate how  the  fit  was  modified  as  the  parameters  change.     This  discussion  is  repeated 
briefly  here. 

From  eqs   (11)  and  (12),    one  can  obtain  a  single  relation  which  needs  to  be  satisfied 
for  a  simultaneous  fit  of  the  experimental  second  virial  and  viscosity  coefficients  at  the 
experimental  temperature,     T.      This  relation  is 

3/2  3/4 

T1(T)  B(T)     =    K  M  B-t  (T-K) 

^3/4  n^^'^^*  ^^*' 

where    K     is  a  constant:     K=137.93xl0  c.g.s.,   units,   and     M     is  the  molecular 

weight.     The  left-hand- side  of  eq  (13)   is  a  function  of  the  experimental  temperature,     T, 
whereas  the  right-hand- side  is  a  function  of  the  reduced  temperature,     T*.     The  equation 
is  used  as  follows:     We  tabulate  the  right-hand- side  as  a  function  of    T*    for  a  chosen 
potential.     At  a  given  experimental  temperature,     T    ,     the  left-hand- side  is  evaluated 
from  experiment.      By  interpolation  in  the  values  tabulated  for  the  right-hand- side,    a  value 
of    T*,      say     T    *,     is  obtained  for  which  the  right-hand- side  of  (13)  equals  the  left-hand- 
side.     Using     e/k  =  T    /T    *    we  thus  find     e/k    at    T       and  also  the  corresponding  value 
of    a     (or    b    )     from  either  eq  (11)  or   (12).     In  a  graphical  sense,    the  pair  of  values     a 
(or    b    )     and     e/k    can  be  represented  by  a  point  in  a  graph  having     a     (or     b    )     and     e/k 
as  axes.      This  point    r  .presents   the   simultaneous   solutLon  of  eqs    (11)  and   (12)  and 
necessarily  corresponds  to  a  simultaneous  fit  of  the  viscosity  and  second  virial  coeffi- 
cients at  the  experimental  temperature     T    .      One  can  determine  points  for  each  of  several 
experimental  temperatures,     T    ,     T    ,     T    ,     T    .     If  the  data  are  exact  and  the  potential 
perfect,     a    and     e/k    values  obtained  at  all  temperatures  are  identical  and  there  results 
a  single  point  in  a  plot  of     e/k    versus     a    as  a  function  of    T.      In  an  actual  calculation 
process,   however,   an  area  is  traced  out  since,    obviously,    the  data  are  not  exact  and  the 
potential  is  not  perfect.      Thus,   both     e/k    and     g    are  tennperature  dependent.     See 
figure  2.     If,   however,    the  potential  is  reasonable,    i.e.,    if  one  could  use  it  to  fit  both 
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second  virial  and  viscosity  coefficients  simultaneously  within  experimental  error  with  the 
same  paranaeters  at  all  temperatures,  then  the  scatter  in  figure  Z  would  be  small  and  of  a 
size  consistent  with  the  scatter  in  the  experimental  data. 

Application  to  Argon 

h 

We   illustrate  how  this  method  works  for  argon.     [We  are  restricted  here  to  the  low 
and  insensitive  temperature  ranges   since  the  high  temperature  range  for  the  second  virial 
coefficient,   as  given  by  eq  (2a),    is  inaccessible  if  one  requires  data  of  the  necessary 
precision  [l  l] .  ] 

The  m- 6  Potential:     The  method  was  first  applied  to  the  thr  ee-paranneter  potentials 
of  the  m-6  family  (corresponding  to  the  value    y  =  0     in  the  four-parameter  nn-6-8  function 
family).     Figure  3   shows  the  curves  for     m  =  9,    12,    17  and  24.     As  already  described,   and 
as  sketched  in  figure  2,    each  of  the  points  on  these  curves  corresponds  to  a   solution  of 
eq   (13)  at  a  given  experimental  temperature.     Notice,    (a)   no  curve    is  reduced  to  a  comnaon 
solution  point  independent  of  temperature;  or,    more  realistically,    there   is  no  small  area 
whose  size  could  be  described  as  a  point  consistent  with  the  experimental  error;   (b)  there 
is,    however,    a   "turn  around"  with  respect  to  temperature  as  a  function  of     m.      For 
example,    examine  the  relative  positions  of  the   solutions  for  the   exper innental  tennperatures 
129.56  K  and  143.  16  K  for  the  9-6  and  the  24-6  functions  [18].     From  this  one  can  surmise 
that,    if  a  common   solution  point  were  to  exist,    it  would  have  to  be  found  for     m    less  than 
24. 

The  m-6-8  Potential:     We   indicated  in  Section  2.  1  and  figure  1   that  increasing  the 
Y     paraineter   (which  represents   the   strength  of  the   inverse  eight  attraction)   effectively 
tends  to  increase     nn.     For   instance,   an  11-6  potential  with    y  =  2.0   is   similar  to  a  12-6 
potential  with    y  =    I'O-      Thus,    based  on  the  curves  for  the  m-6   potential  shown  in  figure 
3,   a  commion  solution  point  cannot  occur  with  the  m-6-8  potential  for     m  >  24     (and, 
possibly,   based  on  an  exannination  of  figure  3,    cannot  occur  for     m  >  17)   since  rerjipving 
the  restriction    y  =  0   is   equivalent,    in  this   sense,    to   increasing     m.      This   initial  crude 
selection  process   is  helpful   since   it  results   in  a  reduction  in  the  number  of  m-6-8 
functions  which  need  to  be  tested  for  argon.     In  keeping  with  this,    we  concentrated  our 
attention  on  potentials  with     m    values  between  9  and  18  and  show  the  particular  results 
obtained  for     m  =    10,    11,    12  and   13  for  various  values  of    y     in  figures  4  and  5.      We  found 
that  a  "turn  around"   is  observed  for  the  12-6-8  and  13-6-8  potentials  at  high  values  of    y. 
This  turn  around  indicates  an  upper  limit  for    y    associated  with  each     m,    i.  e.  ,    for  the 
12-6-8,     y  2  2.0    and  for  the   13-6-8    y   <  1.5.     We  also  found  that  the  common  solution 
area  becanae  progressively  larger  as     m     increased  beyond     m  =   12  no  matter  what  value 
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of    Y     was   selected.     Based  on  this  latter  behavior,    therefore,   we  can  further  restrict  the 

parameter  range  and  say  that,   for  the  m-6-8  potential,     m     must  be  less  than  12  for  argon. 

Closer  examination  showed,    with  very  little  doubt,    that  the  graph  of     e/k    versus    b       had 

o 

the   smallest  common  solution  area  for  the  11-6-8  with    y  =  3.0,    i.e.,    this   is  the 
potential  for  which  the  best  simultaneous  fits  are  obtained  for  the  viscosity  and  second 
virial  coefficients  for  argon  over  the  widest  temperature  range. 


From  within  the  common  solution  area  obtained  in  an     e/k    versus     b       plot  for  the 

3     ° 
11-6-8  potential  with    y  =  3.0,    we  then  chose  the  parameters     b     =45  cm    /mol 

o 

[a  -  3.29Z  A],   and     e/k  =   153  K,    table   1.     We  have  therefore  been  able  to  satisfy  require- 
ment (2)  for  argon  and  with  parameters  which  fit  the  viscosity  data  of  argon 
simultaneously  at  high  and  at  low  temperatures.     This   ls  a  very  strong  argument  in  favor 
of  using  this  potential  for  argon. 

One  might  argue  that,    because  we  now  have  four  parameters,   a  unique  solution  to 
the  problem  of  producing  a  simultaneous  fit  to  the  viscosity  and  second  virial  coefficients 
is  not  possible.      In  other  words,   although  we  find  the   11-6-8,     y  -   "i.O  potential  to  be  the 
best  potential,    we   might  equally  well  have  chosen  another  as  best,    e.g.  ,    the  10-6-8, 
Y  =  4.  5   (say)  or  the  9-6-8,     y  =  ^'      Collision  integrals  were  not  available  to  check  on  this 
but,    in  any  case,    if  multiple  solutions  occur,    they  most  probably  indicate  that  the  data  are 
not  sufficiently  precise  to  make  a  truly  definitive  choice  of  the  potential  possible.     The  data 
are,    however,    sufficiently  precise  to  reduce  the  number  of  possibilities  quite  drastically. 
For  instance,   we  have  already  ruled  out  any  solutions  for     m  >  13,    regardless  of  the  value 
of    Y' 

3.3     Deviation  Curves  for  Argon    _ 

Given  the  potential     m  =   1 1  ,     y  =  3.  0,    with     a  =  3 .  29  A  (1    A  =   1 0~         m)     and 
e/k  =   153  K,    we  constructed  deviation  curves  for  the  dilute  gas  viscosity  coefficients 
[13-16,    19-23]  and  the  second  virial  coefficient  [17],   figures  6  and  7.     We  regard  them  as 
satisfactory. 

4.     Further  Checks  on  the  Use  of  the  11-6-8   (y  =  3.0)  Potential  for  Argon 

Thermal  Conductivity:     The  statistical  mechanical  expression  for  the  thermal 
conductivity  coefficient  of  a  dilute  n:ionatomic  gas,     \'     involves  the   same  collision 
integrals  as  does  the  viscosity  coefficient: 


,  25  V^  mkT  % 

32  2     (2,  2)>.-  m 

■n  o    ii  (T='=) 
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^  /  15      k     ^ 

4      m      ' 

3 

where     c    '/m  =  —  k/m     is  the  translational  specific  heat  per  molecule.     A  check  on  thermal 
V  Z 

conductivity  thus   serves  as  a  cross  check  on  the  viscosity  and  thermal  conductivity  data 
and  not  necessarily  as  a  check  on  the  potential  function.     Figure  8  gives  the  deviation 
curve  for  the  thermial  conductivity.     Data  are  taken  from  references  [Z4-35]. 

Self- Diffusion  Coefficient  D:     This  coefficient  is  of  more  importance  than  the 

thermal  conductivity  as  a  check  on  the  potential  because  a  collision  integral  other  than 

(2    2)* 
Q     '  (T*)     enters  into  the  theoretical  expression.     This  expression  is: 

_  5  TT  mkT  ,,  r. 

P°    =     8  2     (1,1)*  ^  ^''^ 

T^  o    Q     '  (T*) 

where     p     is  the  mass  density  of  the  gas  and     Q     '  (T*)     is  the  collision  integral  for 

diffusion.     Figure  9  gives  a  plot  of  deviations  of  self- diffusion  data,    references  [36-38], 
and  values  calculated  from  eq  (15)  using  the   11-6-8   {y  =  3.0)  potential. 

Isotopic  Thermal  Diffusion  Factor,   a^:     Since  combinations  of  collision  integrals 
are  required  by  the  kinetic  theory  expression  for  this  property,    it  is  especially  sensitive 
to  the  potential  function.     According  to  Kihara  [3], 

a      =    a    '  a  +  6)  (16) 

o  o 


where 


with    6     given  by 


,     _      15(6C»  -   5)   (2A*  +  5) 
^'o       "     2A>1=  (16A*  -    12B*  +  55) 


(17) 


2A*  r  1    /7(5  -  6C*)  +  A*l*  \       /(35/8)  +  28A*  -  6F* 


7:  [-4 


(35/4)  +  7A*  +  4F*     L  2V  5  +  2A*  /       \  21A* 

(18) 
5     r  7     /5  -   6C*  \       31* 


[„*. 


7     L  5     \5  -f  2A*  /        10 
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where 

I*    =     8  -   8E* 
and 

(35/4)  -   3B=:=  -  bC^< 

«*  =  TT-G^^ (19) 

The  potential  function  enters  into  the  equation  for    q'       via  the  collision  integrals  which 
here  are  expressed  as  ratios,  A*,    B-!',    C-~l%    E*  and  F*    given  by 

A*     .     q(2'^)*/q(1'1)% 

B*    =     (5q(1'2)>._^^(1,3)*^/^(1,1)*^ 

C*    =     n^    '    '    /n*    '    '     ;  (20) 

Figure  10  shows  a  plot  of    a       determined  from  eq  (l6)  compared  to  the  data  of  Paul 

at  al.    [39].     Since  thermal  diffusion  data  are  difficult  to  obtain  with  precision,   we  regard 

the  fit  indicated  by  this  figure  as   satisfactory. 

5.     Tests  on  the  Validity  of  the  11-6-8   (y  =  3.  0)  Potential 

Our  results,   as  illustrated  by  the  deviation  curves,   give  us  confidence  that  we 
have  indeed  reached  an  important  goal,    that  is,    we  now  have  a  useful,    relatively  simple, 
analytical  model  potential  function  for  argon.     As  we   shall  show  later,    this   same  model 
potential  is  also  satisfactory  for  the  other  heavy  rare  gases  and  very  reasonable  for 
simple  polyatomic  molecules.     However,    it  is  interesting  to  check  on  how  valid  our  argon 
potential  is;   i,  e.  ,    to  see  how  close  it  connes  to  fulfilling  requirements  other  than  those 
associated  with  the  fitting  of  equilibrium  and  transport  data  [99].     For  example,    the 
repulsive  part  of  the  model  potential  can  be  compared  with  results  from  molecular  beam 
scattering  experiments  and  the  coefficients  of  the  attractive  part,    i.e.  ,     C,     and     C^, 

D  O 

with  independent  semitheoretical  quantum  mechanical  calculations  [101].     Independent 
checks  on  the  values  of     e/k    and     a    are  also  possible.     For  exannple,    the  value  of    o     can 
be  compared  to  that  obtained  from  the  crystal  lattice  [5]  and  the  value  of     e/k    can  be 
compared  with  a  value  found  appropriate  in  a  theory  of  the  liquid  state  [40].     In  both  these 
latter  cases,   however,    one  has  to  introduce  the  concept  of  non-additivity  of  the  pair 
potentials  which  introduces  an  extra  uncertainty  [43], 
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5.  1     Repulsion:     The  Potential  for  r  <  a 

The  repulsive  part  of  the  pair  potential  for  argon  at  small  separations  has  been 
estimated  from  molecular  beam  scattering  data  [41]  and  we  plot  a  composite  curve 
obtained  from  these  data  as  the  heavy  solid  line  in  figure  11.     The  curve  for  the  11-6-8 
(Y  =  3.0),   was  calculated  and  is  shown  as  a  lighter   solid  curve.     The  two  curves  begin  to 
separate  at    r  ~  2.  8  A  or   §  (r)  ~  0.  1   eV.     [Also  shown  as  a  matter  of  interest  are  the 
curves  for  the  12-6  function,   long  used  as  a  good  first  approximation  potential  for  the 
description  of  high  temperature  argon  data,   and  for  the  18-6  which  is  typical  of  all  three 
parameter  potentials  which  can  be  used  to  describe  low  temperature  argon  data  [3,17].] 
Taken  literally  this  comparison  indicates  a  need  for  a  potential  at  small    r     which  is  much 
softer  than  the  11-6-8.     This  conclusion  must  be  taken  with  caution,   however,   because 
repulsive  contributions  do  not  dominate  the  interaction  potential  in  this  region.     For 
example,    our  potential  for  argon  has  a  magnitude  of  approxinnately  1   eV  at  r*  =  0.75.     For 
this  value  of    r*,     the  attractive  inverse  eighth  term  has  a  nnagnitude  which  is  38%  of  that 
of  the  inverse  eleventh  and  the  attractive  inverse  sixth  ternn  has  a  magnitude  which  is  2% 
of  the  inverse  eleventh.     Thus,    if  our  potential  for  argon  is  at  all  valid,    it  follows  that  the 
analysis  of  scattering  data  in  terms  of  pure  repulsion  can  produce  misleading  results. 

An  idea  of  the  error   introduced  by  neglecting  attraction  can  be  estimated  as 
follows.     At  each  value  of    r*,     let  us,    in  the  style  of  the  scattering  experiments,    represent 
our  potential  function  by  a  pure  inverse  power  repulsion.     Thus 

/ 

|(r*)     =    A/r*"^ 

If  this  were  a  true  representation  of  the  function,   an  effective  repulsive  exponent  could 

.        r  -I-    d  f 
then  be  obtained  simply  by  calculating     m     =  — -  - — -  .      But,   for  example,   at    r*  =  0.  75 

'  $    dr* 

this  leads  to  an  effective  value  for  our  potential  of     m'  =   13:     different  fronn  the  value 
m  ::=  11     associated  with  the  repulsive  branch  of  the  complete  potential. 

5.2    Attraction:    A  Comparison  of  the  Dispersion  Coefficients 

The  expressions  for   the  dispersion  coefficients,      C,     and     C    ,     for  the  n:i-6-8 
potential  function  are,   frotn  eq  (10): 

C/     =    --^  [m  -  Y(m-8)]  a     d  (21) 

6  (m-6) 

C_     =     e  y  o     d  (22) 
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For  argon  we  have  found  that    m=ll,     y  =  3,0,      e  =   153Kergsand,     ct  =  3.Z9x10        cm. 

Also,     d  =   1 .  14462.     Using  these  values  in  eqs   (21 )  and   (22)  one  obtains     C ,  s  21   x   10 

ergs/cm       and     C     =2000  x  10  ergs/cna    .     Senaitheoretical  quantum  nnechanical  calcu- 

lations,   however,    indicate  that     C,  =  60  x   10  ergs/cm     and     C     =  120  x   10  ergs/cnn 

6  8 

[42].     The  coefficient    C,     is  believed  to  be  known  quite  accurately  (within  5%),     C.     is 

2  ^ 

known  somewhat  less  well,    but  is  about  twice     C    r    , 

6 

Obviously,    our   11-6-8  potential  for  argon  is  not  giving  the  correct  value  for  the 

inverse  sixth  power  coefficient     C,  .     Let  us  then  see  what  happens  to  the  data  correlation 

if  the     C,     of  our  potential  is  forced  to  have  a  value  approximately  equal  to  60  x  10 

6 
ergs/cnn     .     In  so  doing,    we  effectively  fix  one  of  the  four  parameters.     This  we  choose 

to  do  by  fixing    y    for  a  given     m.     From  the  simultaneous  fits   shown  in  figure  5  we  select 

m  =  11,   and    y  =  2    which  gives     C,   =  55  x   10  ergs/cm     --    sufficiently  close  to  theory. 

Corresponding  best  values  of    a    and     e/k    are  found  from  the  graph  to  be  3.356   A  and 

137  K,    respectively.     Deviation  curves  for  the   second  virial  and  viscosity  using  this 

11-6-8   (y  =  2.0)  potential  are  presented  in  figure   12.      One  sees  that  these  curves  are  not 

unreasonable,    that  is,    the  fit  is  within  three  percent,    except  for  the   second  virial 

coefficient  at  very  low  temperatures.     They  do,   however,    show  systematic  deviations  and 

are  definitely  inferior  to  the  curves  given  in  figures  6  and  7. 

These  results  can  be   sumniarized  in  three  remarks. 

1.  The  numerical  values  of    C,     and     C„     for  the  m-6-8  potential  are  very 

6  o 

sensitive  to  the  choice  of     m    and    y    because  of  the  particular  functional  forms  used, 

i.e.  ,    eq  (21).     The  quantity     a,     which  is  also  part  of  the  expressions,    is  found  from 

experiment.      Since     a     is  raised  to  the  sixth  and  eighth  power  respectively,   uncertainty  in 

it  obviously  leads  to  increased  uncertainty  in     C,     and     C    . 

6  8 

2.  A  convenient  argument  says  that  our  potential  is  not  sufficiently  flexible  to 

fit  data  while  simultaneously  giving  a  good  value  for     C,      (and  to  a  lesser  extent     C    ). 

6  8 

That  our  potential  is  not  flexible  enough  is  obviously  true  --   it  is  only  meant  to  approx- 
imate the  true  potential   (whatever  that  is).     This  argument,    though  necessary,    is  not 
sufficient,   however,    to  recommend  the  need  for  additional  parameters.     We  can  show  this 
by  referring  to  a  recent  potential  function  due  to  Barker  [43].      Because  it  has  additional 
parameters,    Barker's  potential  is  very  flexible  when  compared  with  ours  but,   according 
to  figure  13,    the  attractive  part  of  the   11-6-8   (y  =  2.0)  is  close  to  Barker's  potential  even 

though  his     C       is  fixed  at  the  correct  value.      Consequently,    it  can  be  assumed  that  simply 
6 

adding  parameters  to  the  m-6-8  for  greater  flexibility  might  not  result  in  a     C,     which 
approaches  the  theoretical  value.     In  fact,    it  is  possible  that  the  requirement  of  fitting 


16 


data  as  well  as  we  have  done  in  figures  6  and  7  might  be  incompatible  with  the  requirement 
for  a  value  of    C,      in  reasonable  agreement  with  theory  [99]. 

3.       Since  we  regard  the  data  we  have  fitted  as  reliable,    the     m-6-8    potential  thus 
fails  in  that  it  cannot  give  the  dispersion  coefficients  correctly.     We  must  emphasize, 
however,    that  this  failure  must  not  be  taken  out  of  context.     The  m-6-8  potential,   being 
a  model  function,    can  only  be  approximate.     Its  advantage  is  in  its  utility  and  not 
necessarily  in  its  ability  to  give  correct  theoretical  predictions  for  the  form  of  the  poten- 
tial itself  --  although  it  would  obviously  be  a  very  powerful  property  of  the  potential  if 
this  were  so. 

6.     Application  of  the  m-6-8  Potential  to  Krypton  and  Xenon 

We  have  tested  the  m-6-8  potential  for  krypton  and  xenon  in  the  same  manner  as 
described  for  argon.     For  both  gases,    dilute  gas  viscosity  data  and  second  virial  data 
were  simultaneously  fitted  using  eq  (13).     Plots  were  obtained  similar  to  figures  4  and  5. 

The  data  for  viscosity  were  taken  fronn  references  [14-16,   44]  for  krypton,   and 
references  [14,    15,   45-47]  for  xenon.      The  second  virial  data  for  both  gases  were  taken 
fronn  reference  [17]. 

For  both  krypton  and  xenon  there  was  little  doubt  that  the  use  of  an  m-6-8  with 
m=ll,     Y=3.0     resulted  in  an     e  /k  -  b       plot  with  the  smallest  area,    this  being  our 
criterion  for  the  best  solution  of  eq  (13).     In  this  way  we  obtained  the  best  potential 
parameters  given  in  table  1. 

Viscosity  data  for  both  high  and  low  temperatures  and  second  virial  data  for  all 
temperatures  were  fitted  separately  for  each  gas.     We  confirnned  that  our  choices  of 
paranneters  based  on  the  simultaneous  fits  to  all  data  were  also  valid  for  the  individual 
properties.     Thus  we  found  that  the  nn-6-8  potential  can  be  used  to  satisfy  requirements 
(1)  and   (2)  of  section  1  for  krypton  and  xenon,   just  as  it  could  for  argon. 

Corresponding  States:     We  have,    incidently,    thereby  verified  that,   by  our 
criteria,    the  corresponding   states  principle  is  applicable  to  argon,    krypton  and  xenon. 
The  law  of  corresponding   states  postulates  that  for  fluids  with  naolecules  whose  pair- 
wise  interactions  obey  a  spherical  potential     $(r),     one  can  write 

§  (r)     =     e  f  (r/o) 

where     e    and     a    are  the  energy  and  distance  parameters  respectively,   and    f     is  the 
same  function  for  all  fluids.     Since  the  only  potential  parameters  which  differ  for  each 
of  the  gases  argon,    krypton,   and  xenon  are     a    and     e/k,  corresponding  states      is  valid 
by  this  criterion. 
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The  Coefficient    C,  :     We  have  already  pointed  out  for  argon  that  the  coefficient  of 

D 

the  inverse  sixth  power  term,    when  allowed  to  behave  as  a  parameter  in  fitting  data,   was 
found  to  be  about  1 /2  -    1 /3  of  the  numerical  value  predicted  by  quantum  mechanics.     We 
have  found  this  to  be  also  true  for  krypton  and  xenon:     (theoretical  values  in  brackets  [42]) 
Kr,     C,  ~  43  [~  107]  and  Xe,     C,   ~  100  [~  236],   all  values  times   10       . 

D  D 

6.  1     Deviation  Curves  for  Krypton  and  Xenon 

Figures   14  and  15  give  deviation  curves  for  the  viscosity  and  second  virial 
coefficients  of  krypton  and  xenon,    respectively.     Figures  16  and   17  are  curves  for  the 
thermal  conductivity  coefficients  determined  by  eq  (14)  is  compared  to  experinnent 
[31,35,48,49,50-52].     Figures   18  and  19  are  deviation  curves  for  the  thermal  diffusion 
coefficients  for  the  two  gases  [53,54].     In  figure  20,    we  include  a  deviation  curve  for  the 
recently  measured  [55]   self- diffusion  coefficient  of  krypton. 

7.     Application  to  Polyatomic  Molecules 

The  evidence  presented  in  the  previous  sections  indicates  that  the  mi-6-8  potential 
is  satisfactory  for  the  monatomic  nnolecules.     However,    the  function  would  be  more 
innportant  as  a  practical  tool  if  it  could  also  be  successfully  applied  to  polyatomic 
molecules.     Unfortunately,    this  extension  is  not  straightforward  since  two  very  basic 
complications  are  introduced: 

1.  In  principle,    the  spherically  symmetric  mi-6-8  potential  (or  any  spherical 
potential)  can  never  properly  depict  the  interaction  between  two  polyatonnic  nnolecules; 
the  potential  for  such  molecules  must  be  dependent  on  the  relative  orientation  of  the 
colliding  nnolecules. 

2.  Even  if  one  has  a  potential  which  can  account  for  orientation  dependent 
interactions,    the  statistical  mechanical  expressions  should  allow  that  collisions 
between  polyatonnic  molecules  will  be  inelastic,    that  is,    the  nnolecular  internal  degrees 
of  freedonn  have  to  be  considered. 

The  nnost  common  way  to  resolve  these  problenns   is  to  ignore  them  unless  it  is 
very  obvious  that  one  cannot  do  so  (when  considering  ions,   for  exannple).     This  over- 
sinnplified  approach  can  be  justified  in  a  way  which  we   shall  describe  below.      Let  us, 
therefore,    investigate  the  application  of  the  nn-6-8  potential  to  the  viscosity  coefficient 
and  the  second  virial  coefficient  of  the  polyatomic  gas,    nitrogen,   using  the  same  expres- 
sions and  procedure  as  given  for  the  rare  gases.     The  sources  for  nitrogen  data  are: 
viscosity,    references  [13-16,22,56-58]  and  second  virial,    reference  [17], 


7.  1     Simultaneous  Fit  of  the  Second  Virial  Coefficient  and  Viscosity- 
Coefficient  for  Nitrogen 

We  followed  the  procedure  described  in  section  3.2  by  comparing  the  second  virial 
and  viscosity  coefficients  for  nitrogen  to  calculated  values  based  on  the  m-6-8  using 
eq  (13).     Recall  again  that  we  were  seeking  values  for     m    and    y    which  give  values  for 
b       and     e/k    which  are  essentially  independent  of  temperature  when  allowance  is  made 
for  a  reasonable  amount  of  experimental  scatter. 

Partial  success  was  achieved:     it  was  possible  to  show  that  solutions  for     m  s  13 
did  not  seem  likely,   but  we  were  unable  to  select  that  m-6-8  potential  for     m  ^  12  which 
gave  the  smallest  area  in  a    b     -    e/k    plot.     A  possible  choice  is     m  =   10,     y  =   l-0» 
o  =  3.694  A,      e/k  -  87  K.     Deviation  plots  for  the  viscosity  and  second  virial  coefficients 
determined  fronn  these  paranaeters  are  shown  in  figure  21. 

These  curves,   although  showing  evidence  of  systematic  error,   are  actually  not 
unreasonable,    i.e.  ,   above  200  K  both     T]    and     B    are  fitted  to  within  2%.     We  have, 
nevertheless,   been  unable  to  achieve  results  for  nitrogen  as   satisfactory  as  were  those 
obtained  for  the  rare  gases.     It  thus  appears  that  requirement  (2)  of  section  1   is  not 
satisfied  for  the  polyatonnic  miolecule  nitrogen. 

7.2     Fits  to  Single  Properties:     The  Viscosity  Coefficient 
and  the  Second  Virial  Coefficient 

We  next  follow  the  procedure  outlined  in  section  3.  1,    i.e.  ,    the  correlation  of  the 
viscosity  coefficient  individually.     Here  we  do  have  a  significant  result:    for  nitrogen  as 
for  argon,   a  single  set  of  parameters  was  found  which  enables  the  viscosity  coefficient 
to  be  fitted  simultaneously  at  low  and  at  high  temperatures.     These  parameters  are 
m  =  12,     Y  =  2.0,     a  =  3.54  A,     e/k  =   118  K,    see  table  2  [59].     In  other  words,    the  m-6-8 
potential  can  satisfy  requirenaent  (1)  of  section  3.     Figure  22  gives  the  deviation  curve 
for  the  fit  for  viscosity  and  also  shows  the  corresponding  fit  of  the  second  virial 
coefficient  for  these  same  parameters.     The  latter  plot  indicates  highly  systematic 
deviations. 

The  second  virial  coefficient  can  also  be  fitted  by  itself  and  we  have  done  this  in 
figure  23  using  the  parameters     m  =   12,     y  =  2.0,     a  -  3.60  A,     e/k  =   123  K.     The 
corresponding  viscosity  deviation  curve  is  also  shown  in  figure  23.     Note  the  deviation 
pattern  which  shows  that  the  percentage  deviation  is  essentially  independent  of 
temperature  ! 
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7.3     Fits  for  Methane,    Oxygen  and  Carbon  Dioxide 

We  found  essentially  the   same  results  for  these  gases  as  we  did  for  nitrogen;  an 
attempt  to  fit  both  second  virial  and  viscosity  data  simultaneously  gave  parameters  which 
produced  only  a  moderately  good  representation  of  such  data.     Just  as  for  nitrogen  it  was 
possible  to  fit  the  viscosity  at  both  high  and  at  low  temperature  with  a  single  set  of 
parameters,    see  table  2.      The  parameters  required  for   such  fits  were  found  not  to  be 
suitable  for  the  second  virial  coefficients,    however.     The  converse  was  also  true, 
parameters  obtained  from  a  best  fit  of  the  second  virial  data,    were  not  satisfactory  in  the 
corresponding  fit  for  viscosity.     In  fact,   deviation  curves  were  obtained  similar  to  figure 
23  for  nitrogen. 

Deviation  curves  for  the  viscosity  of  these  gases  are  shown  in  figures  24-26. 
Sources  of  data:  for  oxygen,  see  the  appendix  and  references  [22,47,  56, 64-66],  for 
methane,   references  [56-60,61,63]  and  for  carbon  dioxide,  references  [16,56,62]. 

8.     Kinetic  Theory  of  Polyatomic  Molecules 

The  two  problems  introduced  when  statistical  mechanical  expressions  are  applied 
to  polyatomic  molecules  have  been  discussed  by  Mason  and  his  coworkers  [67].     Let  us 
simply  list  his  conclusions  and  verify  that  our  results  with  nitrogen,    oxygen,    naethane  and 
carbon  dioxide  are  consistent  with  them. 

8.  1     Inelastic  Collisions 

The  effect  of  inelastic  collisions  on  the  statistical  mechanical  expressions  are  as 
follows: 

Viscosity  Coefficient:     The  evidence  indicates  [98]  that  inelastic  collisions 
have  little  effect  on  the  viscosity. 

Diffusion  Coefficient:     Formal  expressions  for  the  diffusion  coefficient  of 
polyatomic  nnolecules  are  available  [98]  but  it  appears  that  they  differ  from  the  corre- 
sponding monatomic  gas  expression  by  a  factor  which  is  very  close  to  unity. 

Thermal  Diffusion:    Approximate  calculations  [91]   indicate  that  the  expression 
for  the  isotopic  thermal  diffusion  factor  as  given  by  eq  (16)   is  not  complete  for  polyatomic 
molecules  and  should  be  modified.     At  the  present  tinne,    it  is  not  clear  how  to  do  this, 
however,    so  that  it  is  common  practice  to  use  eq  (16)  without  alteration. 

Thermial  Conductivity:  The  thermal  conductivity  equation,  eq  (14),  is  definitely 
incomplete  and  has  to  be  modified  for  polyatomic  gases.  We  discuss  this  in  a  separate 
section. 
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8.2    Nonspherical  Potentials 

A  nonspherical  potential  can  be  written  as  the  sum  of  a  spherical  contribution, 
which  is  independent  of  the  relative  orientations  of  the  colliding  molecules,   and  a  non- 
spherical contribution,    which  depends  on  this  relative  orientation.     In  an  exact  calculation, 
the  relative  orientations  must  be  considered  to  change  during  a  collision.     Because  of  this, 
a  proper  determination  of  the  second  virial  coefficient  and  collision  integrals  can  become 
very  complicated.     In  practice,   however,   a  sinnplification  is  made  based  essentially  on 
the  fact  that  a  statistical  averaging  must  eventually  be  performed.     Thus,    one  assumes 
that  the  molecules  collide  with  a  fixed  relative  orientation  and  collision  integrals  are 
calculated  for  each  orientation  independently.     The  final  values  of  the  collision  integrals 
are  then  determined  by  a  statistical  averaging  over  all  these  possible  orientations  [102]. 

A  potential  of  interest  to  us  is  the   12-6  with  a  term  added  to  allow  for  quadrupole- 
quadrupole  interactions: 

$    =     $(12-6)  +   $(quadrupole  -   quadrupole)  (23) 

Collision  integrals  for  eq  (23)  have  been  calculated  [97]  and  were  found  to  be  very  close 
to  the  corresponding  collision  integrals  for  the  12-6  alone.     [However,    the  combinations 
of  the  collision  integrals  required  for  the  thermal  diffusion  factor  are  significantly 
different  from  their  equivalents  for  the  12-6.]     Second  virial  coefficients  for  several  gases 
have  also  been  determined  from  eq  (23).     [Although,    in  this  case,    extra  terms  were  added 
to  account  for  other  nonspherical  effects.]     Unlike  the  collision  integrals,    the  second 
virial  coefficients  were  found  to  be  substantially  different  from  their  counterparts 
calculated  with  the  12-6  alone. 

8.  3     Summary 

Returning  to  our  fitting  procedures:     one  would  conclude  from  the  above  comments 
that  the  viscosity  could  be  fitted  satisfactorily  with  a  spherical  m-6-8  potential  but  a 
simultaneous  fit  of  the  viscosity  and  second  virial  coefficients  then  obtained  would  not  be 
satisfactory.     This  is  what  was  observed.     Work  is  in  progress  to  modify  the  m-6-8 
potential  to  include  nonspherical  effects. 

9.     Thermal  Conductivity  of  Polyatomic  Molecules 

It  is  known  that  the  simple  equation  for  the  thermal  conductivity,    eq  (14),    is  not 
valid  for  polyatomic  molecules.     Several  modifications  have  been  proposed,    but  the  most 
satisfactory  from  a  theoretical  viewpoint  is  the  relatively  simple  expression  obtained  by 
Mason  and  Monchick  who  approximated  the  Wang- Change-de  Boer  theory  for  a  dilute 
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polyatomic  gas  [67].     Mason  and  Monchick  showed  that  the  conductivity  can  effectively  be 
separated  into  two  parts:     a  part  depicting  the  transfer  of  thermal  energy  by  the  transla- 
tional  motion  of  the  molecules,   and  a  part  depicting  the  transfer  of  energy  by  changes   in 
the  internal  energy  of  the  molecules: 

X    =     X'  +  X"  (24) 

where     X'     is  given  by  eq  (14).     Splitting  the  specific  heat  into  a  translational  piece  and 
a  piece  associated  with  internal  degrees  of  freedom 

c       =     c     '  +  c    "  (25) 

V  V  V 

Mason  and  Monchick  write: 


^'     =    -    Tl    (—   -   Aj  (26) 

c" 

^"     =     PD.    ^     (^^  +  a)  (27) 

mt     \    m  / 

were  D.  is  a  diffusion  coefficient  for  internal  energy  and  A  is  a  term  accounting  for 
the  intercharge  of  internal  and  translational  energy.  According  to  Mason  and  Monchick, 
A     is  given  by 

where     Z     is  the  collision  number  for  rotation,    related  to  the  number  of  collisions  required 
for  the  interchange  of  rotational  and  translational  energy. 

If  one  assumes  that    D.         can  be  replaced  by  the  self- diffusion  coefficient  of 

mt 

eq  (13)   (It  is  not  clear  what  difference  this  assunnption  makes)  then  the  formula  for  the 
conductivity  of  a  dilute  polyatomic  gas  becomes 

2m  m  TT\2  T\    /  mZ 

Although  eq  (29)  is  not  new  and  has  often  been  quoted,    it  has  rarely  been  used  to  determine 
numerical  values  of  the  thernnal  conductivity.      The  reasons  for  not  using   it  were  uncer- 
tainty in  the  potential  function  and  the  lack  of  information  for     Z.      Today,    however,    these 
reasons  are  much  weaker:     we  have  presented  an  argunnent  as  to  why  the  nn.-6-8  potential 
can  be  used  to  evaluate  the  viscosity  and  diffusion  coefficients  for  a  simple  polyatomic 
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gas;  and  several  experimental  determinations  of     Z    have  now  been  published.     [Some  of 
these     Z    values  are  based  on  acoustical  measurements  [68]  and  others  on  thermal 
transpiration  studies  [69].] 

We  were  unable  to  find  reliable  values  of     Z    for  carbon  dioxide.     There  are, 
however,   values  of     Z     given  in  reference  [69]  for  nitrogen,    oxygen,   and  nnethane.     Using 
these  and  the  potential  parameters  for  the  viscosity,    table  2,    we  computed  the  thermal 
conductivity  from  eq  (Z9)  for  these  latter  gases.      Corresponding  deviation  curves  are 
given  in  figures  27-29.     Sources  of  data:     for  nitrogen,    references  [27,  28,  32  ,  35  ,  70- 77]; 
for  oxygen,    references  [27,  74,  75 ,  78-85];  and  for   methane,    references  [27,32,34,81,86-88], 
Values  of    c  "     were  taken  from  NBS  Circular  564  [89]  and  reference  [90]. 

Although  there  is  some  evidence  of  systematic  error,    we  are  generally  satisfied 
with  these  curves. 

10.     Thermal  Diffusion  of  Polyatomic  Molecules 

Realizing  that  the  equation  for  the  thermal  diffusion  factor,     a    ,     may  not  be  com- 
plete,  and  that  the  collision  integral  combinations  from  the  m-6-8  potential  nnay  not  be 
adequate  for  polyatomic  molecules,    let  us,    nevertheless,    consider  a  comparison  of    a 
with  experiment  for  nitrogen,    oxygen  and  carbon  dioxide  using  the  unmodified  eq  (16)  and 
the  potential  parameters  obtained  for  viscosity.     The  experimental  results  are  interesting 
because  they  indicate  that    c^       can  be  negative  for  temperatures  below  critical.     Curves 
connparing  theory  with  experiment  are  given  in  figures  30-32.     Data:     nitrogen  [92], 
oxygen  [93],   and  carbon  dioxide  [94].      One  sees  that    a       is  correctly  predicted  to  go 
through  zero  for  oxygen  but  negative  values  are  not  obtained  for  carbon  dioxide  or 
nitrogen.     [See  reference  [95],   however.]     We  do  not  draw  conclusions  from  these  figures: 
In  addition  to  the  theoretical  problems,    the  experimental  data  must  not  be  regarded  as 
reliable  [96]. 

11.      Conclusion 

The  four-paranneter  m-6-8  potential  function  has  been  derived  and  we  have 
dennonstrated  how  it  represents  the  viscosity,    thernnal  conductivity,    second  virial  and 
diffusion  coefficients,   and  the  thermal  diffusion  factor  for  the  gases  argon,   krypton, 
xenon,   nitrogen,    oxygen,    carbon  dioxide  and  methane.     It  appears  the  m-6-8  potential  is 
the  simplest  model  function  which  can  represent  the  above  mentioned  properties  for  the 
rare  gases  with  a  single  set  of  parameters  used  for  each  gas  for  all  properties.     It  can 
also  represent  the  transport  properties  for  the  polyatomic  gases  listed,    but  is  less 
successful  when  simultaneously  applied  to  the  viscosity  and  second  virial  coefficient 
(although  the  second  virial  can  be  fitted  by  itself).     For  the  polyatomic  gases,    it  will  be 
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necessary  to  include  nonspherical  characteristics  of  the  intermolecular  interactions  in  the 
potential.     Work  in  this  direction  is  in  progress.     Work  is  also  in  progress  to  modify  the 
m-6-8  potential  to  allow  for  quantum  effects. 

On  the  basis  of  a  comparison  between  the  repulsive  part  of  the  argon  ra-6-8  (m  =  11) 
potential  and  molecular  beam  experiments,  and  between  values  we  find  for  the  dispersion 
coefficients  and  the  equivalent  values  estimated  from  quantum  mechanics,  it  seemis  that 
the  n:i-6-8  is  too  simiple  to  represent  the  argon-argon  interaction  fromi  a  theoretical 
standpoint.  However,  we  show  that  it  is  not  clear  that  sinnply  adding  additional  parameters 
to  the  m-6-8  will  resolve  this. 
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TABLE  1 
Values  of  the  parameters  of  the  in-6-8  potential  function 
which  fit  both  transport  and  equilibrium  properties 

gas  m  Y  o  e/k 

A  [s  10"  ^°m]  K 

Ar  11  3  3.292  153 

Kr  11  3  3.509  216 

Xe  11  3  3.841  295 


TABLE  2 

Values  of  the  parameters  of  the  m-6-8 

potential  function  which  fit  transport  data 

gas                  m                                             y  a  e/k 

N^                  12                 2  3.54  118 

CH,                 11                3  3.68  168 

4 

0^                  10                 1  3.437  113 

CO^                 14                1  3.68  282 
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Figure  1.   Plots  of  the  m-6-8  potential  computed  from  eq  (10)  to  illustrate  how  the 

attractive  part  varies  with  (a)  an  increase  in  m  at  fixed  y>  (t>)  an  increase 
in  Y  at  fixed  m. 
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Figure  3.   Solutions  of  eq(13)  for  argon  at  several  experimental  temperatures  for  the 

m-6  potential  function.   Note  the  "turn  around','  with  respect  to  temperature, 
of  the  solutions  at  129.56  K  and  143.16  K,  circled  in  the  figure. 
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Figure  10.  Experimental  and  theoretical  [eq(16)]  thermal  diffusion  factors  for  argon. 
Data  from  [39].   We  use  the  same  m-6-8  potential  as  for  the  other  transport 
coefficients . 
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Figure  11.  The  repulsive  potential  at  small  values  of  r.   We  compare  the  11-6-8  (y  =  3.0) 
potential,  found  to  be  satisfactory  for  the  dilute  gas  properties  of  argon, 
with  the  potential  estimated  from  high  energy  scattering  of  beams  of  neutral 
argon  atoms  [41].   Also  shown,  as  a  matter  of  interest,  are  the  curves  for  the 
well-known  12-6  and  18-6  potentials.   [This  latter  represents  the  best  three- 
parameter  potential  fit  for  argon  at  low  temperatures] 
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Figure  12.     Deviation  plots  for  the  second  virial  and  viscosity  coefficients  of  argon 
based  on  the  m-6-8  potential  with  the  inverse  sixth  power  coefficient 
fixed  at  ~  60  X   10-^0  ergs/cm^:     m  =   1 1 ,     y  =  2.  0,     a  =3.356  A, 
e/k  =   137  K.     Data  not  distinguished. 
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Figure  14.   Deviation  curves  for  the  viscosity  and  second  virial  coefficients   of  krypton. 
Data  for  viscosity,  0  [14],  A  [15],  V  [16],   n[44].   Theoretical  values 
determined  for  eqs(ll)  and  (12)  with  the  m-6-8  potential,  parameters  as  given 
in  table  1. 
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Figure  15.   Deviation  curves  for  the  viscosity  and  second  virial  coefficients  of  xenon. 
Data  for  viscosity,  0  [14],  A  [15],  d[45],  V  [46,47].   See  table  1. 
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Figure  18.   The  isotopic  thermal  diffusion  factor  for  krypton  calculated  from  eq(16), 
using  parameters  given  in  table  1,  compared  to  experiment  [53]. 
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Figure  19.   The  isotopic  thermal  diffusion  factor  for  xenon  calculated  from  eq(16),  using 
parameter  given  in  table  1,  compared  to  experiment  [54]. 
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Figure  21.   Deviation  curves  for  the  viscosity  and  second  virial  coefficients  of  nitrogen 

o 
based  on  a  simultaneous  fit.   Parameters:   m  =  10,  y  =   1.0,    a   =   3.694  A  and 

e/k  =  87.0  K.   Data  (not  distinguished):   viscosity  [13-16],  second  virial 

coefficient  [17]. 
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figure  22.   Fits  of  the  viscosity  and  second  virial  coef f f icients  of  nitrogen^based  on  a 
fit  of  the  viscosity  data.   Parameters  m  =  12,  y  =  2.0,  o  =  3.54  A,  e/k  = 
118.0  K.   Viscosity  data,  V  [13],  A  [14],  d[15],  0  [16],  .[22],  O  [56], 
A  [57],  #[58]. 
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Figure  23.   Deviation  curves  for  the  viscosity  and  second  virial  coefficients  of  nitrogen 
based  on  a  fit  of  the  second  virial  coefficient  alone.   Parameters,  m  =  12, 
Y  =  2.0,  o  =  3.60  A,  e/k  =  118.0  K.   Note  (a)  the  unusual  deviation  pattern 
for  the  viscosity,  (b)  that  the  potentials  used  here  and  in  figure  22  differ 
only  by  a  small  amount  in  a. 
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Figure  30.   Isotopic  thermal  diffusion  factor  for  nitrogen  calculated  from  eq(16)  with  the 
potential  parameters  used  for  viscosity.   Data  from  reference  [92]. 
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Figure  31.   Isotopic  thermal  diffusion  factor  for  oxygen  calculated  from  eq(16)  with  the 
viscosity  potential  parameters.   Data  from  reference  [93]. 
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Figure  32.   Isotopic  thermal  diffusion  factor  for  carbon  dioxide  calculated  from  eq(16) 
with  the  viscosity  potential  parameters.   Data  from  reference  [94]. 
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APPENDIX 

Viscosity  of  Oxygen 

Although  oxygen  is  a  conannon  and  important  fluid,    there  are  very  few  reliable  sets 
of  viscosity  data  for  it.     In  fact,    in  our  opinion,    the  reliable  data  lie  in  the  limited  tem- 
perature range  of  about  250  to  400  K.     We  reject  data  outside  the  range  for  the  following 
reasons: 

a)  Above  400  K:    Data  above  400  K  are  available  but  are  due  largely  to  one 
source  -   Trautz  and  co-workers  -   measured  40  years  ago  [64].     But 
Trautz  measured  the  viscosity  for  many  gases  and,   for  gases  other  than 
oxygen  [Ar,   N    ,   He,   air  and  so  on],    modern  experiments  give  results 
which  differ  from  his.      The  results  differ  by  about  1-5%  for  tennperatur es 
up  to  1000  K  --   the  error   increasing  with  temperature.      Since  there  are 
sound  reasons  to  indicate  that  the  more  recent  work  is  nnore  likely  to  be 
correct  for  these  other  gases  [12],   we  have  also  to  suppose  that  the  old 
data  for  oxygen  may  not  be  correct. 

b)  Below  250  K:    Values  reported  below  250  K  were  measured  30  years  ago 
by  Johnston  and  co-workers  [56].     There  appears  to  be  no  more  recent 
work  [100].     But,    the  viscosities  of  Johnston  for  gases  other  than  oxygen, 
differ   slightly  from  more  recent  data.      Since,   here  again,    we  feel  the 
more  recent  data  for  these  gases  to  be  nnore  reliable,    we  have  to  suspect 
the  earlier  data  for  oxygen. 

Discrepancies  between  the  older  and  the  nnore  recent  viscosities  seenn  to  depend 
only  on  the  tennperature  and  not  on  the  gas   (at  least  to  a  first  approxinnation) .     We  checked 
this  by  constructing  a  deviation  curve  between  old  and  nnodern  data  for  argon,   neon, 
nitrogen,    krypton,   xenon,   air,    nnethane  and  heliunn.      The  curve  is   shown  schennatically  as 
figure  Al.      On  the  basis  of  this  figure  we  feel  justified  to  adjust  the  viscosity  results  for 
oxygen  [56,64]  according  to  the  percentage  discrepancy  indicated  at  the  appropriate 
tennperatures.     Adjusted  oxygen  data  are  given  in  table  Al. 

The  adjusted  data,    together  with  data  fronn  references  [22,47,65,66]  were  then 
fitted  to  the  nn-6-8  potential  and  the  result  is  shown  as  a  deviation  curve,   figure  25,    in  the 
nnain  text. 
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TABLE  Al 

Adjusted  Experimental  Viscosities  for  Oxygen 

Temperature  Viscosity 

3 

K  10  g/cm-s 

90.3  0.0679 

118.8  0.0890 
131.3  0.0979 

144.9  0.108 

158.5  0.117 

172.6  0.128 
184.6  0.137 
400.8  0.258 
500.1  0.305 
550.1  0.327 
556.1  0.328 
675.1  0.377 
769.1  0.411 
881.1  0.450 
963.1  0.477 

1102.1  0.521 
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